Skip to content

Abundance speed#34

Open
cgiovanetti wants to merge 2 commits into
mainfrom
abundance_speed
Open

Abundance speed#34
cgiovanetti wants to merge 2 commits into
mainfrom
abundance_speed

Conversation

@cgiovanetti

@cgiovanetti cgiovanetti commented Jun 18, 2026

Copy link
Copy Markdown
Owner

Pull construction of rate_vec out of RHS of ODE solve for abundances to avoid reconstruction of this array at every step in the solve. Now only computes once per call/once per nuclear_rates_q. ~15% speedup of abundance module.

Cara Giovanetti and others added 2 commits June 15, 2026 12:22
Precompute log(T9_vec), log(mu_median_vec), log(expsigma_vec) once at
construction instead of recomputing them on every frwrd_rate_param call.
Algebraically identical (abundances bit-for-bit unchanged); ~5% faster on
60-row rate tables, ~10% on 500-row dense tables (recovers the dense-table
interpolation penalty).

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant